Mott Transition in the Two-Dimensional Hubbard Model 
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Spectral properties of the two-dimensional Hubbard model near the Mott transition are investi- 
gated by using cluster perturbation theory. The Mott transition is characterized by freezing of the 
charge degrees of freedom in a single-particle excitation that leads continuously to the magnetic 
excitation of the Mott insulator. Various anomalous spectral features observed in high-temperature 
superconductors are explained in a unified manner as properties near the Mott transition. 

PACS numbers: 71.30.+h, 71.10.Fd, 74.72.Kf, 79.60.-i 



The metal-insulator transition due to Coulomb interac- 
tions between electrons is called the Mott transition. In 
the insulating phase, the charge excitation has an energy 
gap whereas the spin excitation is usually gapless [l| . On 
the other hand, a metal far away from a Mott insulator 
exhibits free-electron-like behaviors. Thus, the question 
of how the two limits can be reconciled at the Mott transi- 
tion is a fundamental puzzle in condensed-matter physics. 
In particular, in relation to the anomalous properties of 
cuprate high-temperature (high-T c ) superconductors 0- 
4], the Mott transition in a two-dimensional (2D) sys- 
tem has been investigated from various viewpoints [4j-|20| . 
However, its nature remains controversial. 

This Letter shows how the single-particle excitation in 
a 2D metal changes into the spin- wave mode of the Mott 
insulator [l| , through analyses of numerical results on the 
2D Hubbard model according to the exactly known prop- 
erties of a one-dimensional (ID) system 2l|. In addition, 
anomalous spectral features observed in high-T c cuprates, 
such as the pseudogap, Fermi arc, flat band, doping- 
induced states, spinon-like and holon-like branches, as 
well as kink and waterfall in the dispersion relation 0- 
4J , are explained in a unified manner as properties of the 
2D Hubbard model near the Mott transition. The 2D 
Hubbard model is defined by the following Hamiltonian: 
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where Ci a and ni a denote the annihilation and number 
operators of an electron with spin a at site i, respec- 
tively, and means that sites i and j are nearest 
neighbors on a square lattice. The doping concentration 
is denoted by S. We consider the spectral function de- 
fined as A(k,u) = — ilmG(fc, w), where G(k,u) denotes 
the retarded single-particle Green function at momentum 
k and energy uj [7] of the 2D Hubbard model at zero tem- 
perature. We employ cluster perturbation theory (CPT) 
[20} to obtain G{k,ui), using (4 x 4)-site cluster Green 
functions calculated by exact diagonalization. The size 
of clusters treatable with exact diagonalization in CPT is 
larger than that in cellular dynamical mean-field theory 
14|. We consider the properties of hole-doped systems 



Properties in the (0,0)-(tt,tt) direction.— The behavior 
in the (0,0)-(7r,7r) direction [Figs. QJi-l)] resembles that 
of the ID Hubbard model [Fig. Ufa)]. However, the spec- 
tral weight near (0,0) is shifted to lower energies. This 
can be interpreted using the weak-interchain-hopping ap- 
proximation. Up to the first order of interchain hop ping 
t_i_, we obtain G^ 1 (k,u) = G^(k x ,Lj) — tj_(k) [22|. l23j 
using the ID Green function Gm(k x ,u) and t±(k) = 
— 2t± cos k y . In this approximation (called RPA), the 
spectral weight of a chain shifts to higher [lower] energies 



for t ± (k) > [t±(k) < 0] [23L l24j|. Thus, the spectral 
weight near (0,0) [t±(k) m —2t < 0] shifts to lower ener- 
gies, whereas that near (tt/2,w/2) [t±(k) 0] is almost 
unaffected by interchain hopping [Fig. Efb)]. As a result, 
the dispersion relation exhibits a kink between (0,0) and 
(7r/2,7r/2) [Figs. QJi-1)]. This kink can be identified as 
the high-energy giant kink observed in high-T c cuprates 
[H, Hq . In high-T c cuprates, a nearly vertical dispersion 
relation known as a waterfall [§] has been observed just 
below the kink, involving smaller spectral weights than 
those above the kink. The RPA is too simple to explain 
the reduction in spectral weight. Nevertheless, the results 
obtained using CPT [Figs. H{i-l,p)], where higher-order 
t^-corrections are involved within the clusters, show that 
waterfall behavior appears in the 2D Hubbard model [3] . 

The RPA argument above allows us to trace the origins 
of the dominant modes in Figs. [TJ^i— 1) back to those of ID 



with t > for < k y < k x < it without loss of generality. 



2lll25|,l26|. In Figs. [Tfi— k), the mode for u> > originates 
from the upper edge of the spinon-antiholon continuum 
[Fig. Ufa) , dashed red curve] . The mode for — 1 < uj/t < 
and < k x s» k y < it/2 primarily originates from the 
spinon mode [Fig. Efa), dashed-dotted blue curve]. The 
mode for —4 < u/t % —1 primarily originates from the 
holon mode [Fig. Eta), solid purple curve]. The mode 
almost extending over the Brillouin zone for — 6 < u/t < 
—4 originates from the shadow band [Fig. [2fa), solid 
light blue curve]. In Fig. HJl), the modes originating 
from the ID spinon, holon, and shadow-band modes can 
be identified similarly. The mode originating from the ID 
spinon mode and that from the ID holon mode have been 
observed in high-T c cuprates [3]. The mode originating 
from the ID shadow band should also be observed if other 
bands outside the 2D Hubbard model do not intervene. 
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FIG. 1: Results on the 2D Hubbard model for U/t = 8 obtained using cluster perturbation theory, (a-d) A(k,u)t and (e-h) 
A(k)[= J divA(k,oj)] for 5 = 0.2, 0.07, 0.03, and (from above). In (e-h), contributions from to < 0, the lower Hubbard band 
(LHB) for uj > 0, and the upper Hubbard band are denoted by solid blue curves with hatches, dotted red curves, and dashed 
pink curves, respectively. Total spectral weights (solid brown lines) satisfy A(k) = 1 within numerical accuracy, (i-1) Close-ups 
of (a-d) for the LHB. Rightmost panels show A(ui)t[= J dkA(k, ui)t / (2ir) 2 ]. In (1), dotted green line indicates energy of the top 
of the LHB. Bandwidth of the mode originating from the ID spinon mode, pseudogap defined by the flat mode at (7r,0), and 
that by the peak in A(u>) at 5 — are denoted by £(o,o), e (ir,o) an d e P [pink arrows at (0,0), (it,0), and in the rightmost panel], 
respectively, (m-o) A(k,uj)t at uj as for 8 = 0.2, 0.07, and 0.03 (from above), (p) Close-up of (k) near cj = in the (0,0)-(7r,7r) 
direction. Dashed blue curve indicates e2D{k,k)/t = —y/2(v2D/t)cosk, where spin-wave velocity of the 2D Heisenberg model 
V2D ~ 1.18V2J Hi} (J = U 2 /U). Solid green lines show cj = 0. Gaussian broadening is used with standard deviation a = O.lt. 



A characteristic feature of the Mott transition is the 
doping-induced spectral-weight transfer from the upper 
Hubbard band (UHB) to the lower Hubbard band (LHB) 
04E EH- As shown in Figs. QJi-k), the spectral 

weight transferred to the LHB for uj > is primarily car- 
ried by the mode originating from the upper edge of the 
spinon-antiholon continuum. The energy of the mode at 
(7r,7r) [e(7r, 7r)] does not reach zero even in the 5 — > limit 
[Fig. He)], and its spectral weight gradually disappears 
as 5 — > [Fig. Hpi)]: the mode remains dispersing with 
the spectral weight fading away as 6 —> [Fig. [TKp)] . 
This feature is the same as that in ID [Figs. [2jd,h)] [21] 
and contrasts with that of a doped band insulator and 
with that of a Fermi liquid as the effective mass m* — > oo 
where the effective bandwidth e*(k) oc 1/m* 0]. 



We consider its relationship to the spin-wave mode at 
5 = Q. As shown in Fig. [ljp), the dispersion relations 
of the mode for uj > and the mode originating from 
the ID spinon mode for uj < in the (0,0)-(7r,7r) direc- 
tion behave as €2D(k,k) — —^/2v2D cos k near the Mott 
transition in the large U/t regime, where V2D denotes the 
spin- wave velocity of the 2D Heisenberg model 27(. In 
fact, as shown in Fig. EJe), the bandwidth of the mode 
for uj > in the LHB in the 6 — > limit [er-n-^)] an d 
that of the mode originating from the ID spinon mode 
at 5 = [e(o,o)] [Fig- \H})] reasonably agree not only with 
each other but also with v^^d in the large U jt regime. 
Thus, the modes lead continuously to the spin- wave mode 
at 6 = whose dispersion relation in the (0,0)-(7r,7r) di- 
rection is expressed as E(k, k) — |e2n(fe — n/2, k — ir/2)\ 
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FIG. 2: (a) ,4(fc, w)t in the lower Hubbard band (LHB) for 
U/t = 8 at 5 ~ 0.07 in the ID Hubbard model obtained using 
the dynamical density-matrix renormalization group method 
[2lj . Curves indicate dominant modes obtained using exact 
solutions (b) RPA results for t± = t obtained using 

the data in (a). In (a,b), straight solid green lines indicate 
uj — 0. Gaussian broadening is used with standard deviation 
a = O.li. (c) e(n, ir) [energy of the mode for ui > in the LHB 
at (7r,7r)] for U/t = 8 [Figs. [D;i-k,p)] (blue circles). Solid red 
curve indicates a fit. Dashed green curve shows results for 
U — 0. (d) Same as (c) but for ID. Solid blue curve shows 
e(7r)/i for U/t = 8 [21| . In (d,f), e n , and eo are defined in 
the ID Hubbard model similarly to e(7r,7r), and £(o,o)i 

respectively, (e) t^.^/t [e(n,n)/t extrapolated to S — > 0] 
(red circles) and t(o.o)/i [Fig. [Hi)] (blue diamonds). Dotted 
green line indicates \/2«2D /t, where spin-wave velocity of the 
2D Heisenberg model v 2D « 1.18-^/2./ [13 (J = 4t 2 /U). Inset 
shows J/t dependence, (f) Same as (e) but for ID. Solid blue 
curve shows e w /t(= to/t) [2ll |. Dotted green line indicates 
vid /t, where spin-wave velocity of the ID Heisenberg model 
vid — ivJ/2 29]. (g) e p /t (blue diamonds) and t^, o)/t (red 
circles) [Fig. 0^1)] . Dotted green line indicates a fit in the 
large U/t regime, assuming e p , £(^,o) oc J. Inset shows J/t 
dependence, (h) Spectral weight A for uj > in the LHB. 
Blue squares with solid line show 2D results for U/t = 8. 
Red circles with dashed line denote ID results for U/t — 8 
taken from ref. [2l]]. Dotted green line indicates results for 
t = [l(J. Inset shows the Mott gap A between the LHB and 
the upper Hubbard band at S — 0. Blue diamonds show 2D 
results. Dotted red curve denotes ID results [33| . 



[1|, [27| . The continuous evolution to the spin- wave mode 
is consistent with the scaling behavior of spin correla- 
tions @, [13 18 1 . In the extremely large U/t regime, 
the ferromagnetic fluctuation arising from Nagaoka fer- 
romagnetism [28j could dominate the antiferromagnetic 
fluctuation near the Mott transition, an outcome that is 
beyond the scope of the present study. Note that the re- 
lationship to the spin-wave mode is essentially the same 
as that in ID [2l[: in ID, the dispersion relations of the 
mode for w > in the LHB and the spinon mode reduce 
to ei£i(fc) = — vid cos k in the S — > limit in the large U/t 
limit [Figs. [Ha,f)] [2l|, where v\d denotes the spin- wave 
velocity of the ID Heisenberg model [2{|. The modes in 
ID lead continuously to the spin excitations at 8 = 
whose dominant mode shows E{k) = |eip(fc — 7r/2)| [29j |. 

We next discuss hole pocket behavior. In ID, the re- 
gion between the two gapless points at k = w(l — S)/2 and 
tt(1 + 36)/ 2 [intersections of solid purple and pink curves 
with u = in Fig. HJa)] can be regarded as a hole pocket 
2 1| . Its signature should persist near the Mott transi- 
tion even in 2D because t±(k) ~ near (7r/2,7r/2). The 
small intensity at the intersection of the upper edge of 
the continuum [bending back near (7r/2,7r/2)] with to = 
for k x w k y > tt/2 [Figs. [TJj,k,n-p)] might be regarded 
as a signature for hole pocket behavior as in ID. 

Properties near (ir,0).— As shown in Figs. the 
dominant mode near (jr,0) near the Mott transition is 
located below w = 0, whose dispersion relation is anoma- 
lously flat. This anomalously flat dispersion relation has 
been found in numerical simulations 12, lj| 19| and in 
high-T c cuprates Below, we discuss its relationship 
to the pseudogap and Fermi arc. 

Near the Mott transition, the main peak of A(ui)[= 
J dkA{k,oj)/(2n) 2 ] is located below oj = [Fig. 
rightmost panel] . Namely, the spectral weight decreases 
from the peak value as ui — > 0. This feature is called a 
pseudogap. To discuss the pseudogap quantitatively, we 
define two energies: e p as the energy difference between 
the peak in A(u>) and the top of the LHB at 6 = [Fig. 
HJl)], and e( Wl o) as the energy difference between the flat 
mode at (n,0) and the top of the LHB at S = [Fig. 
[Hi)] . Figure [^g) shows that the pseudogap defined in 
A{oj) is primarily due to the flat mode at (7r,0) and that 
the pseudogap is induced by interactions even without 
further-neighbor hopping. In addition, this figure indi- 
cates that the pseudogap is proportional to J (= At 2 /U) 
[e p rj £(7r,o) ~ 0.8 J] in the large U/t regime, implying 
that it will be related to the antiferromagnetic fluctuation 



L2j,|l5|. The pseudogaps defined by A(u) and by the flat 
mode decrease as S increases [Figs, m^k)], because the 
chemical potential is lowered with the flat mode almost 
unchanged. The pseudogap closes at a 5 value where the 
peak in A(u>) or the flat mode crosses w = 0. 

The flat mode also causes Fermi arc behavior in the 
pseudogap regime. In the (0,0)-(vr,0) direction, there 
is no dominant mode crossing lo — in the pseudogap 
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regime, because the flat mode is located below ui = 
[Figs. [TJj ,k)l . Thus, there is no Fermi surface along this 
direction |12Hl5l [. This feature can also be confirmed 
from n(fe)[= |° m duA{k, uj)] [Figs. Hff-h), solid blue 
curves with hatches]: n(k) along this direction in the 
pseudogap regime [Figs. QJf,g), rightmost panels] almost 
saturates to that at 5 = [Fig. QJh), rightmost panel]. 
For the (7r,0)-(7r,7r) direction, two possibilities have been 
discussed: the absence of the mode crossing u> = 0, which 
leads to hole pockets [12Hl4j. and the presence of the 
mode forming the large Fermi surface centered around 
(7r,7r) l~5j . Figures [IJk,o) show intermediate behavior. 
Namely, the mode crossing w = in the (7r,0)-(7r,7r) di- 
rection can be identified as that originating from the up- 
per edge of the spinon-antiholon continuum, whose spec- 
tral weight fades away as S — > 0. As a result, the spec- 
tral weight almost disappears along this direction around 
uj = 0. This leads to the behavior that can be regarded 
as a Fermi arc (disconnected portion of the large Fermi 
surface) 2] [Fig. QJo)]. This behavior contrasts with that 
of the large Fermi surface in the large doping regime [Fig. 
[ljm)]. Fully quantitative explanation for the Fermi arc 
observed in high-T c cuprates requires further study. Also, 
properties in a very small energy scale inaccessible with 
(4 x 4)-cluster CPT are beyond the scope of this Letter. 

Mott gap.— In a (doped) Mott insulator, the Mott gap 
is due to Coulomb interactions, whereas the band gap 
in a (doped) band insulator is due to a chemical poten- 
tial difference between sites or orbitals within a unit cell. 
A deeper insight into Mott physics can be obtained by 
tracing the origin back to ID. In ID, the quasiparticle 
responsible for the Mott gap has been identified as that 
defined by the fc-A string in the Bethe ansatz and called 
the doublon [2l| , which can be regarded as a pair of elec- 
trons ■ By the presence of the doublon, the UHB can 
be distinguished from the LHB 21J. The doublon is not 
a double occupancy because the latter exist in the LHB 
as well as the UHB for U/t < oo. Because the Mott gap 
in 2D can change into that in ID by reducing the in- 
terchain hopping, the quasiparticle that determines the 
energy scale of the UHB in 2D (doublon in 2D) should be 
a descendant of the ID doublon, i.e., a pair of electrons 
rather than a double occupancy. The similarities in the 
spectral weight transferred by doping and those in the 
Mott gap at S = between ID and 2D [Fig. fflh.)] seem 
to support the traceability. Apart from the presence of 
the doublon, we can identify the dominant modes in the 
UHB [Figs. UKa-d)] as in the LHB for u> < 0, noting the 



21|. 



particle-hole symmetry at 5 = 0, as in ID 

Noting that the high-energy magnetic excitations in a 
magnetic field in quasi-lD Heisenberg antiferromagnets, 
which can be regarded as repulsively interacting hard- 
core boson systems, have been explained using two-string 
solutions [24j, |3l|, |32j, we can generalize the concept of 
Mott physics: the generalized doublon (a pair of single 
particles definable by a string with a length of two in ID 



2ll l3ll |32[ or its descendant in higher dimensions [24| ) 
induced by repulsive interactions will be responsible for 
the high-energy states without multi-site (multi-orbital) 
unit cells, in contrast to band theory. The Mott physics 
could also be explored in cold atomic systems. 

Summary.— Anomalous spectral features observed in 
high-T c cuprates, such as the pseudogap, Fermi arc, fiat 
band, doping-induced states, spinon-like and holon-likc 
branches, as well as kink and waterfall in the dispersion 
relation, were explained in a unified manner as properties 
of the 2D Hubbard model near the Mott transition. The 
properties near (tt,0) are characterized by the flat mode 
unlike ID features. The physics of the Mott gap was ex- 
amined by tracing the origin back to ID. The Mott tran- 
sition is characterized by a dispersing mode that leads 
continuously to the spin-wave mode of the Mott insula- 
tor with the spectral weight fading away toward the Mott 
transition due to charge freezing. The loss of charge char- 
acter from the dispersing mode will be a general feature 
of Mott transitions, which contrasts with the transition 
to a band insulator without spin-charge separation. 
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